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I.  INTRODUCTION 


To  create  a  synthetic  aperture  radar  (SAR)  image  using  ground-penetrating  radar 
(GPR)  data,  one  must  know  the  round-trip  distance  between  the  antenna  aperture  and  the 
point  to  be  focused  on.  Given  a  beginning  point  O  in  media  II,  a  “buried”  point  P  within 
media  I,  and  the  relationship  between  the  points,  the  index  of  refraction  n  and  the  total 
“luminal”  distance  R  =  n\R\  +  112R2  must  be  derived  (see  Figure  l).1  If  we  do  not  worry 
about  refraction,  the  distance  may  be  easily  calculated  using  simple  geometry.  When 
refraction  is  taken  into  account,  the  problem  becomes  more  difficult.  When  the  refraction 
is  governed  by  a  lossy  material,  the  problem  becomes  even  more  intractable  and  an 
appeal  to  numerical  methods  becomes  necessary. 


1  The  total  one-way  travel  time,  T  ,  for  a  ray  of  light  is  scaled  by  the  speed  of  light  in  each  respective 
media  T  =  {R\/v\)  +  (Rijvi).  The  physical  situation  may  be  interpreted  as  a  scaling  of  the  respective 
path  lengths,  R  =  n\R\  +  n2  R2,  such  that  T  =  R/c. 
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A  natural  question  is:  Why  would  we  be  interested  in  accounting  for  refraction? 
There  are  two  immediate  advantages  to  doing  so.  First,  correcting  for  refraction  would 
deposit  the  energy  received  by  the  antenna  in  the  correct  image  pixel  (i.e.,  at  the  correct 
physical  distance),  thereby  placing  any  image  anomaly  (i.e.,  a  potential  target  of  interest) 
at  the  correct  depth.  Second,  not  accounting  for  refraction  would  effectively  blur  the 
image  and  degrade  image  fidelity.  The  disadvantage  is  increased  computational 
complexity  and  processing  time.  However,  efficient  software  coding  and  carefully  chosen 
numerical  techniques  can  reduce  the  complexity  and  processing  to  an  insignificant 
percentage  of  the  total  image-processing  algorithm. 

Deriving  the  equations  that  describe  the  optical  path  length  between  the 
originating  point  O  and  the  terminating  point  P  is  best  done  through  multiple 
transformations.  The  purpose  of  these  transformations  is  to  express  the  equations  in  a 
form  easily  understood  with  minimal  mathematical  complexity.  The  price  to  be  paid  is 
measured  in  proliferation  of  equations.  However,  at  the  end  of  the  process,  the  resulting 
equations  are  few  in  number,  comparatively  simple,  and  written  in  terms  of  typical  media 
parameters  (dielectric  constants,  relative  permeabilities,  etc).  Finally,  Sections  II. A  and 
II. B  closely  follow  the  derivations  presented  in  Stratton’s  Electromagnetic  Theory  [1]  but 
with  updated  notation  and  use  of  magnitudes  and  phases  instead  of  separation  into  real 
and  imaginary  parts. 
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II.  SNELL’S  LAW 


A.  COMPLEX  ANGLES  AND  SNELL’S  LAW 


Given  wave  numbers  k2  e  R  and  k ,  e  C ,  define  k\  =  \k\\el%.  The  standard  fonn  of 
Snell’s  Law,  A:isin#i  =  fesin#?,  will  result  in  an  imaginary  0\. 

k-, 


sin  #,  = 


K 


sin#., 


cos#,  =±  1 1 - e  2‘x  sin2  #, 


(1) 

(2) 


with  the  sign  chosen  by  the  boundary  conditions  and  the  requirement  that  the  field  is 
finite  at  all  points  in  space.  This  is  a  problem  because  physically,  the  refraction  angle 
cannot  be  complex.  Moreover,  we  require  a  real  geometric  angle  to  focus  SAR  imagery. 
To  facilitate  our  understanding  of  the  relationship  between  the  complex  angle  0\  and  the 
physical  refraction  angle  we  expand  cos  0\  using  complex  notation  such  that 

cos#,  =  pew  (3) 


with  a  solution  found  for  p  and  i//by  squaring  and  using  (2).  The  solution  to  the  above  is 
straightforward 

k  2 

p1  cos  2  y/  =  1 - T_sin2  #2  cos2j 


k  ^ 

p1  sin2^  =  1 - ^ysin2  #2  sin2/ 


(4) 

(5) 


and  it  is  trivial  to  show 


P  = 


/  \  1/4 

k  2  Jr  4  ^ 

1  -  2  sin2  #2  cos  2 %  +  -^-sin4 


v  & 


ur  =  —  arctan 
2 


K 


L7sin2#2sin2^ 
kx  “  -  k2  sin2  #2  cos  2^ 


(6) 


(V) 
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B.  ELECTROMAGNETIC  WAVES  IN  BOTH  MEDIA 


We  assume  in  media  II  an  electromagnetic  plane  wave  is  present  such  that  it 
intersects  the  boundary  between  lossless  media  II  and  lossy  media  I  (%  ^  0)  at  an  incident 
angle  02.  Let  Ej,  Er,  and  Et  be  the  magnitude  of  the  incident,  reflected,  and  transmitted 
electromagnetic  fields,  respectively.  The  total  electromagnetic  wave  present  will  be 

II:  Eje^-^  +  Ege^'7-^ 

I:  E Te^2'r~wi] 

in  areas  I  and  II.  An  electromagnetic  wave  propagates  normal  to  the  planes  of  constant 
phase,  k  -r  =  constant.  For  all  angles  of  incidence  except  &  =  0,  the  planes  of  constant 
phase  will  not  necessarily  be  parallel  to  the  planes  of  constant  amplitude.  Consequently, 
the  field  within  medium  I  will  not  constitute  a  standard  plane  electromagnetic  wave.  For 
simplicity,  align  ki  =  k]  (sin  Opt  +  cos$,z)  such  that 

k  ■rl=  | kx  \e'x  (sin  0xx  +  cos  9xz)  ■  {xx  +  zz)  then  using  ( 1 ) 

=  z^L^pe1^^  +  xk2smd2  .  (8) 

=  i\k}  \psm{x  +  y/)z  +  \kt  \pcos(%  +  i//)  +  k2  sin6*2x 

Define 

p{02 )  =  p{02  )sin(/  +  if/{d2 )) 
q{d2 )  =  p{02  )cos(z  +  y/{92 )) 
so  the  transmitted  electromagnetic  wave  will  be 

e~\k\  | p-  |?z+^2  sin  e2x-coi) 

in  lossy  media  I.  Constant  phase  planes  in  media  I  will  be  defined  by 

z  cos  "%+  x  sin  E,  =  constant 

from  which  the  true  form  of  Snell’s  Law  is  revealed  for  lossy  media 

.  „  k-i  sin  02 

srn  q  -  .  ~  -1  . . -  . 

2J\kl  “  q{02  )2  +  k\  sin2  02 


(9) 

(TO) 

(ID 

(12) 

(13) 


We  will  find  it  profitable  to  explicitly  note  the  above  equation’s  “cousin”: 

kxq02 

cos#  =  . 

J  kx  I"  q{02  f  +  k2  sin2  02 


Whereas  the  phase  velocity  of  light  in  the  lossless  media  is 
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(15) 


a  careful  look  at  (11)  will  reveal  the  phase  velocity  of  light  in  the  lossy  media  in  our 
geometry  of  interest  (i.e.,  two  semi-infinite  regions  separated  by  a  plane)  is  dependent  on 
the  angle  of  incidence, 

f  (16) 

-Jj&j  |“  q{02  )2  +  k 2  sin2  02 

We  can  define  the  ratio  of  the  velocities  as  the  ratio  of  the  real  index  of  refraction  nx  {02 ) 
in  media  I  to  the  index  of  refraction  «2  in  media  II 

n[02  v2  sin#,  \\kx  |2  f  V  .  2 

-^  =  ^  =  —^=  L^rq(02)  +  sin  02  .  (17) 

n2  vx  sin£  k2 

After  some  manipulation  we  arrive  at  the  following  summarized  results: 


12  \  2 

p(A)=  i-^ysin2^  +4^_sin2^ 


^)4aretanLVSm^S:n25 


2  Ik,  -  k22  sin2  022% 


q{02 )  =  p{&2  )cos(z  +  y(o2 )) 


tan  e  _  ^2  sin  (92 


C.  SANITY  CHECK:  REDUCING  TO  STANDARD  FORM 


Let  x  =  0,  then  (17)-(21)  reduce  to 


y/{02 )  =  ^-arctan(o)  =  0 


Pfe)  = 


k2  I2  1  i - 

1  — |-sin2  02  =  — tJ/c?  -  k\  sin2  02 

k\  k] 


q(&i)  =  y^ki  -k2sm202 
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(25) 


nx  (02)  _  k\  _  nx 

Yl 2  k2  Tl*) 

Given  the  above  results  and  (13)  we  immediately  see  that 

.  „  G  sin  ft 

sin£  =  , 

\\K  f  <l{@2  )2  +  ^2  sin2  02 

=^>  sin£  =  —sin 02 
fh 


(26) 

(27) 


which  is  the  standard  Snell’s  Law  for  refraction  at  the  boundary  between  two  lossless 
media. 


D.  RELATING  WAVE  NUMBERS  TO  TYPICAL  PHYSICAL  PARAMETERS 


After  reducing  ( 1 7)— (2 1)  to  the  standard  form  of  Snell’s  law  when  %  =  0,  we 
relate  k\  and  /o  to  typical  known  parameters.  First,  assume  medium  II  is  air  so  that  ki  ~ 
oWUb/4))  =  cole.  Second,  there  are  a  number  of  standard  descriptions  for  a  complex  wave 
number.  If  medium  I  is  completely  described  by  a  complex  permittivity  e ,  and  a  real 
permeability  //i ,  then 


kx  =  c o-,Jsxjux  =  a>$fce- 

|^i  |  —  Mi 


(28) 


tan  x  =  tan —y . 


(29) 


In  general,  the  permittivity,  penneability,  and  conductivity  of  medium  I  are  all 
functions  of  frequency.  When  exact  solutions  are  required,  the  appropriate  s\ (co),  /u\(co), 
and/or  <j\ (co)  must  be  inserted  into  (28)— (29).  In  practice,  however,  for  many  materials  ju\ 
is  very  near  juo,  and  s\  and  <7\  are  only  weakly  dependent  on  frequency. 


The  physical  characteristics  of  many  lossy  materials  are  often  described  in  the 
literature  through  the  use  of  the  relative  (real)  penneability  //,  =  (jU\/ju 0),  the  real  dielectric 
constant  K  =  (Re{<s)/£b),  and  the  loss  angle  8  such  that  tan  8  =  (Im{£)  /Re  {£•}). 
Therefore,  the  previous  equations  for  p{ d2)  and  i/X  /%)  become 


p(Pi)  = 


1- 


sin2  0, 


km, 


+  sin  8 


r  J 


1/4 


(30) 
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i//(82)  =  ^-arctan 


sin"  6 2  sin  28 


2 k/jr  -  2  sin  02  cos"  8 


(31) 


A  little  algebra2  will  provide 


cos 


\s+v(e,) 

V2 


1 

-  +  ■ 
2 


1 


1  + 


tan"  8 


( ^  sin2  02  ^ 


v 


km, 


r  J 


(32) 


which,  combined  with  the  description  of  the  lossy  media  (i.e.,  media  I)  using  //r,  K,  and  8, 
describes  the  entire  physical  phenomena  of  refraction  with  the  three  equations: 


q{o 2)  = 


^  sin2  02  ' 


v 


km, 


+  sin2  8 


) 


1/4 


-|l/2 


1 

— +  - 
2 


1 


1  + 


tan2  8 


f  sin2  ^ 

1  - 


v 


Km, 


J 


e  _  sin#.,  cos# 


(d2)_  r kMr 


cos  8 


q{ej  +  sin1 02 . 


(33) 


(34) 


(35) 


McLaurin  expansion  of  (32)  will  give  some  physical  insight  as  to  the  strength  of  lossy 
refraction  compared  with  non-lossy  refraction. 


cos 


\s+v(e,) 

V2 


.  1  tan  8  3  tan4  8 

1 - ^  + 


,  sin"  02 

1  - 


KMr 


2  16  (  sin2  6  34 


1- 

v  km,. 


1/2 


(36) 


2  See  Appendix  A  for  a  short  derivation. 


7 


where  we  immediately  see  eos(c5>/2  +  yAOi))  =  1  and  sin//  =  [sin 6h/ ( VAT //,)]  for  8 «  tv! A. 
Therefore,  materials  with  small  loss  angles  can  be  treated  as  essentially  lossless  and  the 
angle  of  refraction  calculated  from  the  familiar  form  of  Snell’s  Law. 

E.  ILLUSTRATION  BY  EXAMPLE 

For  a  soil  mixture  of  50%  sand,  35%  silt  and  15%  clay  with  a  volumetric  moisture 
content,  mv,  of  0.05,  experimental  measurements  give  K  «  5  and  tan  S~  1/10  over  the  300 
MHz-1.3  GHz  frequency  range.  Increasing  the  moisture  content  to  mv  =  0.20  results  in  K 
«  18  and  tan  8  ~  2/18  over  the  same  frequency  range  [2],  In  Figure  2  we  see  that  the 
effect  of  a  nonzero  loss  angle  is  extremely  minor.  For  all  values  of  the  angle  of  incidence 
02,  the  refraction  angle  is  approximated  extremely  well  by  the  standard  form  of  Snell’s 
law.  At  the  limit/%  =  n/2,  the  difference  between  the  case  K  =  18,  tan  8=  2/18  and  K  = 
18,  tan  8=  0  is  less  than  4  x  10  4  radians.  In  Figure  3  we  see  the  lossy  material  has  an 
indices  of  refraction  ratio  weakly  dependent  on  the  angle  of  incidence. 


Figure  2.  Snell’s  Law:  Lossless  versus  Lossy  Refraction 
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Figure  3.  /n2  :  Lossless  versus  Lossy  Refraction 


Therefore,  as  seen  in  Section  II.D,  materials  with  small  loss  angles,  8  «  ti/4,  can 
be  treated  as  essentially  lossless.  The  ability  to  treat  weakly  lossy  materials  as  lossless 
will  greatly  reduce  the  computational  complexity  of  focusing  GPR  and  simplify  the 
relevant  equations  enough  to  allow  for  a  closed-form  solution. 


Figure  4  illustrates  an  effect  that  will  be  of  great  interest  to  us  in  the  following 
section,  namely  the  limiting  refraction  angle  28,  as  the  angle  of  incidence  approaches  the 
grazing  angle: 


lim  28,  — >  2  arctan 


cos  8 


q(o 2)  V 


r  j 


(37) 


This  effect  will  compete  with  the  reduction  in  the  wavelength,  A,  =  A{)/ n\  ,  to  determine 

the  theoretical  resolution  of  the  GPR  system.  The  figure  plotted  assumes  the  simplest 
case  of  ju\  =  //o  and  frequency-independent  K  and  8.  For  small  loss  angles,  the  effect  of 
absorption  on  the  limiting  angle  of  refraction  is  practically  nonexistent.  It  is  only  for  loss 
angles  greater  than  approximately  id  A  that  the  effect  becomes  noticeable.  Evidently,  the 
effect  of  the  dielectric  constant  is  significantly  more  important  in  determining  the  limiting 
refraction  angle. 


9 


Figure  4.  Refraction  Limit  2^ for  <%  =  ?d2 
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III.  GPR  BACKPROJECTION  AND  SNELL’S  LAW 


A.  THEORETICAL  SITUATION 


Figure  5.  Theoretical  Situation 

Focused  GPR  backprojection  requires  knowledge  of  the  round-trip  time, 
t  =  (nt  Rt  +  n2R2)/c  between  the  points  in  the  antenna  aperture  and  the  points  in  space  we 

will  focus  on.  We  begin  by  assuming  that  we  know  both  the  coordinates  of  the  antenna 
aperture  O  and  the  point  P  we  wish  to  focus  on:  (jc2,  yi,  zo)  and  (xi,  y  1,  z\).  We  will  further 
assume  we  have  a  right-handed  coordinate  frame  where  positive  (increasing)  values  of  z 
point  into  media  I.  For  media  II,  we  take  the  standard  case  of  112  =  1  (air),  and  initially 
assume  media  I  is  a  lossy  soil  mixture  of  50%  sand,  35%  silt,  and  15%  clay  (K  =  5,  tan  8 
=  1/10).  However,  because  of  the  insignificant  loss  angle  of  many  soils  and  the 
significantly  increased  computational  complexity  of  treating  the  soil  as  lossy,  we  will 
eventually  treat  the  soil  as  lossless  and  reap  significant  savings  in  processing  time. 

From  Figure  5,  standard  geometry  gives: 
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H  =  -z2 

(38) 

D  =  Z] 

(39) 

=  V(x2  ~  X1  )2  +  (y2  -  yl  )2 

(40) 

r2=  h 

cos  6 2 

(41) 

R,=  D 
cos^ 

«.  L-a  L-H  tan^ 

tanq  = - = - - 

D  D 

(42) 

.cosfi'  L-H  tan6*0 

q(o2)i  A>,  D 


from  which  we  define  (and  require) 


F(62,K,S,iurD,L,H)  =  DS^k_  +H  tan  ft  -1  =  0 


#2)  V  K/U  r 


(43) 


with  the  distance  n2i?2  +  /?,  /?,  given  by 


d  +  D  11  Dnl 
n2R2  ■+■  /7|  R^  — - h 


H  D  Icosd 


cos  6 2  cos^  cos  6 2  q{02 )  \  K/j, 


^-q(0j  +sm2  02 
cos  o 


(44) 


Our  task  is  to  find  the  roots  of  (43)  that  lie  within  the  space  (h  e  [0,  n!2\,  namely  the 
angle  of  incidence  given  the  coordinates  of  the  origin  and  the  point  we  wish  to  focus  on. 
This  will  allow  us  to  calculate  our  true  interest,  n\  R]  +n2R2.  The  nature  of  q((h)  will 

make  finding  a  closed-form  solution  difficult  for  lossy  materials  but,  for  lossless 
materials  ( 8  =  0),  the  solution  to  (43)  can  be  written  in  closed  fonn  using  the  quartic 
formula  to  find  the  roots  of  the  fourth-order  polynomial: 

H2  (l  -  Kjur  )x4  -  2HL(\  -  Kjur  )x2  +  [d2  -  H2Kjur  +  L2  (l  -  Kjur  )Jx2 
+  2  HLKjurx-L2Kjur 

where  62  =  arctan  x.  In  solving  (45),  care  must  be  taken  in  choosing  the  correct  root,  x  e 
[0,  00].  Since  many  soils  are  characterized  by  small  loss  angles,  equation  (45)  can  often  be 
used  in  place  of  the  significantly  more  complicated  (43),  resulting  in  decreased 
computational  complexity.3 


3  See  Appendix  B  for  further  discussion. 
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B.  NUMERICAL  ANALYSIS 


Assume  that  the  synthesized  aperture  will  focus  at  the  point  (0,  0,  3)  and  the 
antenna  moves  along  a  grid  5  cm  above  the  surface.  Using  Matlab  R14  we  are  able  to 
solve  for  the  angle  of  incidence  given  the  aperture  and  “target”  coordinates,  and  the 
results  are  illustrated  in  Figure  6  for  our  test  soil  mixture.  The  details  of  the  figure  match 
physical  intuition:  It  is  symmetric  about  the  origin  and  disappears  at  the  origin. 


Figure  6.  Angle  of  Incidence  Given  P  at  (0,  0,  3)  and  O  at  (x,  y,  -5) 


To  qualitatively  understand  how  important  refraction  is  in  determining  the 
luminal  path  length  R  =  n2Rl  +n2R2,  we  will  measure  the  difference  between  the  true 

path  traveled  and  the  geometric,  straight-line  path  defined  by  n\r\  +  n^r 2. 


/7,/i  +  n2r2 


+  H 


(46) 


The  maximum  one-way  difference  between  the  luminal  path  length  and  the  geometric 
path  length  is  approximately  8  cm  (or  0.27  ns)  when  the  antenna  phase  center  is  at  (-25, 
0,  -5)  cm  and  is  focusing  at  the  point  (0,  0,  7)  cm.  It  is  difficult  to  intuitively  decide  how 
much  of  an  effect  these  differences  would  have  on  GPR  data  processing.  On  one  hand, 
between  500  MHz-6  GHz  an  8  cm  difference  will  be  approximately  0.13  to  1.6 
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wavelengths,  which  can  introduce  serious  phase  errors  when  summing  the  signal  over  the 
aperture.  On  the  other  hand,  when  the  antenna  is  operated  perpendicular  and  near  the 
surface,  most  of  the  difference  appears  only  at  points  where  most  of  the  signal  received 
by  the  antenna  would  probably  be  of  small  amplitude  compared  with  signals  received 
from  points  with  no  significant  need  of  refraction  correction.  However,  antennas  with 
main  lobes  of  large  angular  extent  or  geometries  where  the  antenna  is  oblique  to  the 
surface  may  increase  the  relative  signal  amplitudes.  Therefore,  the  relative  importance  of 
correctly  accounting  for  refraction  can  only  be  determined  on  a  per  antenna  basis  where 
Figure  7  is  weighted  by  the  antenna  beam  pattern  and  due  consideration  is  given  to  the 
geometry  of  the  physical  system. 


Figure  7.  One-Way  Difference  with  P  at  (0,  0,  d)  and  O  at  (x,  0,  -5) 
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IV.  CONCLUSION 


In  summary,  we  have  derived  the  Snell’s  Law  equations  describing  the  refractory 
behavior  of  electromagnetic  waves  intersecting  a  boundary  between  lossless  media  and 
lossy  media  under  the  assumption  of  the  existence  of  electromagnetic  plane  waves  in 
media  II.  Furthermore,  using  Snell’s  Law,  we  have  explicitly  derived  the  equations 
needed  to  describe  the  geometry  of  a  standard  GPR  system-target  package  and 
quantitatively  calculated  the  expected  correction  for  a  test  case  utilizing  realistic  soil 
parameters.  However,  the  validity  and  urgency  of  correctly  accounting  for  the  refractive 
behavior  is  strongly  dependent  on  the  parameters  of  the  GPR  antenna  utilized  and  the 
relative  orientation  of  the  antenna  to  the  boundary. 
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APPENDIX 


A.  DERIVATION  OF  cos 


\s+v(e,) 


Consider  the  square  of  the  term  we  are  interested  in  and  expand 

cos2  \S  +  ^  +  ^-[cosc7  cos2^  -  sinc7  sin2^] . 


Then  from  equation  (3 1)  we  immediately  have 

_  sin2Jsin26,T 

tan  2^  = 


A 


2 KjU,.  -  2  cos"  8  sin  02  B 


(47) 


(48) 


and  therefore 


cos  2y/  = 


sin2^  = 


B 

Ja2+b2 

A 

V A2 +B 2  ' 


Simple  algebra  gives 

A2  +B2  =4cos 2  8 
from  which  we  may  easily  see 

cos^cos2^  = 


(Kn,  -  sin"  6 2 )”  +  ( K/ur  )2  tan°  8 


K/jr  -  cos1  8  sin2  d2 


sin  £sin2^  = 


-  sin2  d2)  +  ( K/urf  tan2  8 

sin2  8  sin2  02 

(A//,  -  sin2  6 2 )  +  (Kjur  )2  tan2  8 


Therefore, 


cos" 


-S+v 


1  1 

=  -  + 


Kjur  -  sin1  02 


^  ^  y] ( Kju r  -  sin2  6 2  y  +  (K/j,.  f  tan2  8 


(49) 

(50) 

(51) 

(52) 

(53) 

(54) 
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B.  QUARTIC  POLYNOMIAL  DESCRIPTION  FOR  LOSSLESS  MATERIALS 

If  8=  0,  then  (43)  reduces  to 

D  Sin#,  1  ttl  n  r  n  /c 

— Y — ,  +  H  tan  02  -  L  =  0  (5i 

cl\@i)  jKfir 

£>sin<9,  tt.  n  ,  a  ,r> 

,  -  +  i/tanff,  -  L  =  0 .  (5 

^Kjur-  sin2# 

Squaring  both  sides  and  replacing  tenns  gives 

°  Sm  °2  +  (H tan ff,  -  Z)2  +  2(7/ tan 02-L)  ,  Dsm^-  =  p  (5; 
Kjur  -sm~02  ^  -  JK/u,.  -  sin2  02 


D  sin“  02  t\2  T x  Dsm62 


KjUr  -  sin2  6*2 


D2  sin2  6*2  -  (//  tan  02-L )2  (/<//,  -  sin2  6*2 )  =  0  . 

It  will  be  advantageous  to  make  the  substitution 

tan  #7  =  x 

2 

— >  sin2  #  =  — * — . 

x2 +  1 

Then  after  expansion  and  term  collection,  we  have 

D2x2  -  (Hx  -L)2  [Kjur  (x2  +  l)  -  x2  ]  =  0 

and  we  can  easily  derive  the  quartic  polynomial 

H2(  1  -  Kjur  )x4  -  2HL{\  -  Kjur  )x3  +  [d2  -  H2Kjur  +  L2  (l  -  Kjur  )]x2 
+  2 HLKjurx  -  L2Kjur 


Practical  use  of  (63)  will  require  care  taken  in  choosing  the  correct  root,  x  e  [0, 00) . 
Although  not  explored  in  this  paper,  finite  group  theory  may  shed  some  light  in 
understanding  the  behavior  of  the  roots  of  (63)  as  a  function  of  D,  H,  L,  K,  and  pir. 
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